#Install packages
install.packages("haven")
trying URL 'https://cran.rstudio.com/bin/macosx/big-sur-arm64/contrib/4.3/haven_2.5.3.tgz'
Content type 'application/x-gzip' length 1132879 bytes (1.1 MB)
==================================================
downloaded 1.1 MB
The downloaded binary packages are in
/var/folders/8t/cf05ypcs56z1d0qh3cl43l880000gn/T//RtmpgK7aj1/downloaded_packages
#Load packages
library(data.table)
library(tidyr)
library(haven)
#Add data WV6
WV6_data <- load("/Users/laurabazzigher/Documents/GitHub/risk_wvs/data/dataset/WV6_dataset_wave_5_6/WV6_Data_R_v20201117.rdata")
WV6_data <- WV6_Data_R_v20201117
print(WV6_data)
#Name the variables
data6 = WV6_Data_R_v20201117[, c("V1", "V2", "V10", "V23", "V70", "V71","V72","V73", "V74B", "V75", "V76", "V77", "V78", "V79", "V147", "V149","V198", "V199", "V201", "V208","V209", "V210", "V240", "V242", "V248", "V229", "V57", "V58")]
names(data6) = c("v1_wave", "v2_country", "v10_happyness", "v23_satisfaction", "v70_creativity", "v71_money", "v72_security", "v73_goodtime", "v74b_help", "v75_success", "v76_risk", "v77_proper", "v78_environment", "v79_tradition", "v147_religion", "v149_hell", "v198_gov_benefits", "v199_fare_public_transportation", "v201_cheating_taxes", "v208_beating_wife", "v209_beating_children", "v210_violence", "v240_sex", "v242_age", "v248_education", "v229_employed", "v57_married", "v58_children")
#Exlusion of participants which no info about sex, age, education,
employment status, marital status, or number of children
data6 = subset(data6, v76_risk > 0 & v240_sex > 0 & v242_age >0 & v248_education > 0 & v229_employed > 0 & v57_married > 0 & v58_children >= 0)
Scatterplot of Age vs. Risk with ggplot
data_df <- data.frame(Age = data6$v242_age, Risk = data6$v76_risk) # Convert data into a data frame.
ggplot(data_df, aes(x = Age, y = Risk)) +
geom_point(color = "blue", size = 1) +
labs(title = "Scatterplot of Age vs. Risk",
x = "Age",
y = "Risk")

Histogramm des Alters
hist(data6$v242_age,
main="Häufigkeitsverteilung des Alters",
xlab="Alter",
ylab="Häufigkeit",
col="blue")
#Number of countries
length(unique(data6$v2_country)) # what is wrong?
#Number of participants
nrow(data6)
#Range of age
range(data6$v242_age)
#Table about gender
table(data6$sex) #is this correct???
#Categorical age variable (for example, to plot response frequencies
by category)
data6$agecat <- NA # Create an empty column 'agecat'
data6$agecat[data6$age < 20] = "15-19"
data6$agecat[data6$age >= 20 & data6$age < 30] = "20-29"
data6$agecat[data6$age >= 30 & data6$age < 40] = "30-39"
data6$agecat[data6$age >= 40 & data6$age < 50] = "40-49"
data6$agecat[data6$age >= 50 & data6$age < 60] = "50-59"
data6$agecat[data6$age >= 60 & data6$age < 70] = "60-69"
data6$agecat[data6$age >= 70 & data6$age < 80] = "70-79"
data6$agecat[data6$age >= 80] = "80+"
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKI0luc3RhbGwgcGFja2FnZXMgCmBgYHtyfQppbnN0YWxsLnBhY2thZ2VzKCJoYXZlbiIpCmBgYAoKI0xvYWQgcGFja2FnZXMgCmBgYHtyfQpsaWJyYXJ5KGRhdGEudGFibGUpCmxpYnJhcnkodGlkeXIpCmxpYnJhcnkoaGF2ZW4pCmxpYnJhcnkoZ2dwbG90MikKYGBgCgoKI0FkZCBkYXRhIFdWNgpgYGB7cn0KV1Y2X2RhdGEgPC0gbG9hZCgiL1VzZXJzL2xhdXJhYmF6emlnaGVyL0RvY3VtZW50cy9HaXRIdWIvcmlza193dnMvZGF0YS9kYXRhc2V0L1dWNl9kYXRhc2V0X3dhdmVfNV82L1dWNl9EYXRhX1JfdjIwMjAxMTE3LnJkYXRhIikgCldWNl9kYXRhIDwtIFdWNl9EYXRhX1JfdjIwMjAxMTE3IApwcmludChXVjZfZGF0YSkKYGBgCiNOYW1lIHRoZSB2YXJpYWJsZXMgCmBgYHtyfQpkYXRhNiA9IFdWNl9EYXRhX1JfdjIwMjAxMTE3WywgYygiVjEiLCAiVjIiLCAiVjEwIiwgIlYyMyIsICJWNzAiLCAiVjcxIiwiVjcyIiwiVjczIiwgIlY3NEIiLCAiVjc1IiwgIlY3NiIsICJWNzciLCAiVjc4IiwgICJWNzkiLCAiVjE0NyIsICJWMTQ5IiwiVjE5OCIsICJWMTk5IiwgIlYyMDEiLCAiVjIwOCIsIlYyMDkiLCAiVjIxMCIsICJWMjQwIiwgIlYyNDIiLCAiVjI0OCIsICJWMjI5IiwgIlY1NyIsICJWNTgiKV0KCm5hbWVzKGRhdGE2KSA9IGMoInYxX3dhdmUiLCAidjJfY291bnRyeSIsICJ2MTBfaGFwcHluZXNzIiwgInYyM19zYXRpc2ZhY3Rpb24iLCAidjcwX2NyZWF0aXZpdHkiLCAidjcxX21vbmV5IiwgInY3Ml9zZWN1cml0eSIsICJ2NzNfZ29vZHRpbWUiLCAidjc0Yl9oZWxwIiwgInY3NV9zdWNjZXNzIiwgInY3Nl9yaXNrIiwgInY3N19wcm9wZXIiLCAidjc4X2Vudmlyb25tZW50IiwgInY3OV90cmFkaXRpb24iLCAidjE0N19yZWxpZ2lvbiIsICJ2MTQ5X2hlbGwiLCAidjE5OF9nb3ZfYmVuZWZpdHMiLCAidjE5OV9mYXJlX3B1YmxpY190cmFuc3BvcnRhdGlvbiIsICJ2MjAxX2NoZWF0aW5nX3RheGVzIiwgInYyMDhfYmVhdGluZ193aWZlIiwgInYyMDlfYmVhdGluZ19jaGlsZHJlbiIsICJ2MjEwX3Zpb2xlbmNlIiwgInYyNDBfc2V4IiwgInYyNDJfYWdlIiwgInYyNDhfZWR1Y2F0aW9uIiwgInYyMjlfZW1wbG95ZWQiLCAidjU3X21hcnJpZWQiLCAidjU4X2NoaWxkcmVuIikKYGBgCgoKI0V4bHVzaW9uIG9mIHBhcnRpY2lwYW50cyB3aGljaCBubyBpbmZvIGFib3V0IHNleCwgYWdlLCBlZHVjYXRpb24sIGVtcGxveW1lbnQgc3RhdHVzLCBtYXJpdGFsIHN0YXR1cywgb3IgbnVtYmVyIG9mIGNoaWxkcmVuCmBgYHtyfQpkYXRhNiA9IHN1YnNldChkYXRhNiwgdjc2X3Jpc2sgPiAwICYgdjI0MF9zZXggPiAwICYgdjI0Ml9hZ2UgPjAgJiB2MjQ4X2VkdWNhdGlvbiA+IDAgJiB2MjI5X2VtcGxveWVkID4gMCAmIHY1N19tYXJyaWVkID4gMCAmIHY1OF9jaGlsZHJlbiA+PSAwKQpgYGAKCiMgU2NhdHRlcnBsb3Qgb2YgQWdlIHZzLiBSaXNrIHdpdGggZ2dwbG90CmBgYHtyfQpkYXRhX2RmIDwtIGRhdGEuZnJhbWUoQWdlID0gZGF0YTYkdjI0Ml9hZ2UsIFJpc2sgPSBkYXRhNiR2NzZfcmlzaykgIyBDb252ZXJ0IGRhdGEgaW50byBhIGRhdGEgZnJhbWUuCgpnZ3Bsb3QoZGF0YV9kZiwgYWVzKHggPSBBZ2UsIHkgPSBSaXNrKSkgKwogIGdlb21fcG9pbnQoY29sb3IgPSAiYmx1ZSIsIHNpemUgPSAxKSArCiAgbGFicyh0aXRsZSA9ICJTY2F0dGVycGxvdCBvZiBBZ2UgdnMuIFJpc2siLAogICAgICAgeCA9ICJBZ2UiLAogICAgICAgeSA9ICJSaXNrIikKYGBgCgojIEhpc3RvZ3JhbW0gZGVzIEFsdGVycwpgYGB7cn0KaGlzdChkYXRhNiR2MjQyX2FnZSwgCiAgICAgbWFpbj0iSMOkdWZpZ2tlaXRzdmVydGVpbHVuZyBkZXMgQWx0ZXJzIiwKICAgICB4bGFiPSJBbHRlciIsCiAgICAgeWxhYj0iSMOkdWZpZ2tlaXQiLAogICAgIGNvbD0iYmx1ZSIpCmBgYAoKCiNOdW1iZXIgb2YgY291bnRyaWVzCmBgYHtyfQpsZW5ndGgodW5pcXVlKGRhdGE2JHYyX2NvdW50cnkpKSAjIHdoYXQgaXMgd3Jvbmc/CmBgYAoKI051bWJlciBvZiBwYXJ0aWNpcGFudHMKYGBge3J9Cm5yb3coZGF0YTYpIApgYGAKCiNSYW5nZSBvZiBhZ2UKYGBge3J9CnJhbmdlKGRhdGE2JHYyNDJfYWdlKSAKYGBgCgojVGFibGUgYWJvdXQgZ2VuZGVyCmBgYHtyfQp0YWJsZShkYXRhNiRzZXgpICNpcyB0aGlzIGNvcnJlY3Q/Pz8KYGBgCgoKCiNDYXRlZ29yaWNhbCBhZ2UgdmFyaWFibGUgKGZvciBleGFtcGxlLCB0byBwbG90IHJlc3BvbnNlIGZyZXF1ZW5jaWVzIGJ5IGNhdGVnb3J5KQpgYGB7cn0KZGF0YTYkYWdlY2F0IDwtIE5BICAjIENyZWF0ZSBhbiBlbXB0eSBjb2x1bW4gJ2FnZWNhdCcKCmRhdGE2JGFnZWNhdFtkYXRhNiRhZ2UgPCAyMF0gPSAiMTUtMTkiCmRhdGE2JGFnZWNhdFtkYXRhNiRhZ2UgPj0gMjAgJiBkYXRhNiRhZ2UgPCAzMF0gPSAiMjAtMjkiCmRhdGE2JGFnZWNhdFtkYXRhNiRhZ2UgPj0gMzAgJiBkYXRhNiRhZ2UgPCA0MF0gPSAiMzAtMzkiCmRhdGE2JGFnZWNhdFtkYXRhNiRhZ2UgPj0gNDAgJiBkYXRhNiRhZ2UgPCA1MF0gPSAiNDAtNDkiCmRhdGE2JGFnZWNhdFtkYXRhNiRhZ2UgPj0gNTAgJiBkYXRhNiRhZ2UgPCA2MF0gPSAiNTAtNTkiCmRhdGE2JGFnZWNhdFtkYXRhNiRhZ2UgPj0gNjAgJiBkYXRhNiRhZ2UgPCA3MF0gPSAiNjAtNjkiCmRhdGE2JGFnZWNhdFtkYXRhNiRhZ2UgPj0gNzAgJiBkYXRhNiRhZ2UgPCA4MF0gPSAiNzAtNzkiCmRhdGE2JGFnZWNhdFtkYXRhNiRhZ2UgPj0gODBdID0gIjgwKyIKYGBgCgoKCgoKCgoKCgoK